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Abstract. We present the results of a monitoring campaign of the double quasar SBS 1520+530 at Maidanak observatory from 
April 2003 to August 2004. We obtained light curves in V and R filters that show small-amplitude Am as 0.1 mag intrinsic 
variations of the quasar o n time scales of about 100 days. The data set is consistent with the previously determined time delay 
of At = (130 ±3) days bv lBurud et all (20021. We find that the time delay corrected magnitude difference between the quasar 
images is now larger by (0.14 ± 0.03) mag than during the observations bv lBurud et alJ 1 120021) . This confirms the presence of 
gravitational microlensing variations in this system. 
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1. Introduction 

The broad absorption line (BA L) quasar SBS1520+530 
(z q = 1.855) was discovered by IChavushvan et alJ dl997l) 
as a gravitationally lensed double quasar with an angu- 
lar sej)aratioji_ojM^56 i _TT^ was detected 
bv lCrampton. Schechter & Beuzid ill 99 8) using infrared adap- 
tive optics i maging at the Canada-France-Hawaii Telescope. 
Faur e et all (|2002) obser ved the lensing g alaxy with the 
Hubble Space Telescope. iBurud et all £2002) (in the follow- 
ing B02) finally succeeded in obtaining both the redshift 
Zgai = 0-717 (con sistent with absorption lines first found by 
Cha vushvan et alJ Il997l) of the lensing galaxy with a Keck 
obervatory spectrum and the time delay of At = (130 ± 3) 
days between the two quasar images using monitoring data 
from the Nordic Optical Telescope (NOT). For this, B02 ob- 
tained an almost gapless lightcurve of the object of about 
800 days between February 1999 and May 2001. An almost 
continuous coverage of the light curve of SBS 1520+530 is 
made possible by its high declination. Further photometry ob- 
tained at Maidanak obs ervatory on the syste m was published 
bv lZheleznvak. Sergeev & Burkhonov! d2003h . 

Interestingly, B02 found that by simply shifting the light 
curve of image B backward by 130 days and correcting for the 
magnitude difference Am = 0.69 mag of the images made the 
quasar light curves align only approximately. They obtained 
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a better match by allowing for an additional linear trend, and 
by correcting for f aster variations using an iterative scheme 
jBurud et all 1200 lb. B02 interpret these additional variations 
to be probably due to microlensing variability. 

SBS 1520+530 thus is one of the prime targets for mi- 
crolensing studies since it provides at once the prospects for 
long, uninterrupted light curves with a known, relatively short 
time delay and known microlensing variations. For this reason 
we decided to continue the optical monitoring of this system. 



2. Observations 

We report here on our observations of SBS 1520+530 (Fig.Q) 
with the 1 .5 m AZT-22 te lescope in Maidanak, Uzbekistan 
( Ehgamberdiev et al., 2000), between April 25 and October 12 
in 2003, and again between January 16 and August 22 in 2004. 
In these intervals we observed the object almost daily whenever 
the weather permitted it. The V and /^-filters provide photom- 
etry in the Johnson-Cousins system. We observed on average 
four frames per night with an exposure of 3.5 min in V and 3 
min in R with a median seeing of 1.07 arcsec. 

In total we obtained 116 nights in the V band and 191 nights 
in the R band. After selecting for good seeing (better than 1.4 
arcsec full-width at half maximum) and low sky-background 
we finally used 80 nights in the V-band and 123 nights in the 
fi-band. The observations were made with the BROCAM CCD 
detector with a pixel scale of 0.26 arcsec. Bias correction and 
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Fig. 2. fi-band zoom of the central part of SBS 1520+530. The 
field size is 30 arcsecx25 arcsec. North is up and East is to the 
left. The bright foreground star next to the two quasar images A 
and B is saturated in the centre. SI and S2 are also fo reground 
stars. The lens galaxy (m« ga i as 21.6 mag. ICrampton et alJ 
1998) is too faint to be seen in this image. 



Fig.l. fl-band image of SBS1520+530 obtained on May 7, 
2003. The quasar images, the reference stars S3 and S4, and 
the star we use for the PSF are labelled. The field size is 2.8 
arcminx3.5 arcmin. The area marked with the box is shown in 
Fig. El 

fiat fielding of the images were done using standard IRAF soft- 
ware. 

3. Photometry 

A special property of the SBS 1520+530 system is its location 
within 14 arcseconds of a bright 12th magnitude star (Fig. [2}. 
This is great for adaptive optics studi es of the system bec ause 
a bright reference star is at hand JCrampton et all fl998t> . For 
photometry, however, the bright halo and the diffraction pattern 
caused by the star on the CCD needs to subtracted carefully. 

In order to do this, we followed the method described by 
IZheleznvak etail §003); we extracted the western half of the 
star and subtracted it from the eastern part where SBS 1520 is 
situated. This procedure also efficiently subtracts the light due 
to the horizontal diffraction spike that extends towards the dou- 
ble quasar. 

Photometry on the quasar components A and B was per- 
formed using the DAOPHOT package dStetsoni Il987h . We 
chose this method because it is well-suited to the mildly asym- 
metric and time- variable AZT-22 point spread function (PSF). 
Because of the variable PSF i t was not po s sible to use the imag e 
subtraction technique jAlard & Luptonl Il998l lAlardl l2000l) . 
The quoted error bars are the standard lcr errors determined 
by DAOPHOT. Since there are four point sources (the quasar 



components A and B, and the stars SI and S2) crowded in 
a rather small region, we had to fit the positions and magni- 
tudes of all four components at the same time. Absolute quasar 
magnitudes were calibrated using the brightness of the refer- 
ence star S3 (see Fig.Q) in the V-band (my = (17.37 ± 0.02) 
mag) and j?-band {rriR = (17.18 + 0.02) mag) determined by 
IZheleznvak et alJ J2003h . The star we used as a template to 
model the point spread function in DAOPHOT is marked with 
"PSF" in Fig. [U 

In our analysis we ignore the presence of the lensing galaxy 
near image B becau se it is too faint to be de tected in our images 
(jw«, g ai = 21.6 mag. lCrampton et aflll99 81 If all of the galaxy 
light contributes to our magnitude estimate of quasar image B, 
the measured brightness would increase by a constant offset of 
0.08 mag (assuming m^p = 18.8 mag). 

4. Results 

4.1. Light curves 

We present the results of the R-bmd photometry of 
SBS 1520+530 in Fig. □ and Table l 1 . The light curves of 
the two quasar components A and B are plotted together with 
the light curve of an additional reference star S4 (see Fig. \l\ 
m.R = 18.16 mag). This shows that both S3 (which is used for 
the absolute magnitude calibration) and S4 do not vary. 

Both quasar components show low-amplitude Am a; 0.1 
mag variations on time scales of about 100 days. As we will 
show in the next section, the overall similarity of the two quasar 



1 Table 1 is being made available electronically at CDS. Columns 
1-4 contain the A and B magnitudes with errors, Column 5 contains 
the Julian Date of the observations. 
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Fig. 3. /?-band light curves of the two quasar images A and B in 
SBS 1520+530 and the reference star S4. For clarity the mag- 
nitude for image B was shifted by -0.4 mag, S4 was shifted by 
+2.4 mag. The dotted vertical line indicates the day with the 
Julian Date-2452000=880 when the mirror was cleaned. The 
dashed horizontal line shows the magnitude of the reference 
star S4 11154 = 18.16 mag. 
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Fig. 4. V - R colour curve for the quasar images A (top) and B 
(bottom). The best-fitting value for each quasar image is plotted 
with a dashed line. 



image light curves is due to a gradual brightness decrease of 
the lensed quasar. On the day with Julian Date-2452000=880 
the telescope mirror was cleaned. This led to a large sensitivity 
improvement that visibly improved the accuracy of the quasar 
brightness measurements, especially for the fainter image B. 



In Fig.0]and Table 2 2 we show the difference between our 
V-band and fl-band light curves of SBS 1520+530. We find an 
average V - R colour m v - m« = 0.15 mag for image A and 
my -m R =0.18 mag for image B. We do not find any evidence 
for significant colour variations during our observing interval. 
The small difference A(V - R) « 0.03 mag of the V - R colour 
between the quasar images indicates the presence of a small 
level of differential reddening along the light paths. 

4.2. Time delay 

A wide variety of algorithms have been developed for the de- 
termination o f time delays in gravi t ational l ens systems (e.g., 
i Kundic et all Il997t iBurud et all l200lt Idil-Merino et all 
2002]). We choose here a strightfoward linear interpolation 
scheme because our light curve mainly consists of two fre- 
quently sampled observing intervals. The unfrequently sam- 
pled gap in the middle cannot be confidently interpolated with 
any method. In detail we use the following recipe: 

1 . The B light curve is shifted by the time delay Af to be tested. 

2. One of the two light curves is linearly interpolated to match 
the observing dates of the other light curve. 

3. Only gaps that are less than 40 days (including the gaps 
between Julian Date-2452000=1000 and 1100) are inter- 
polated, no difference is calculated for larger gaps because 
they can introduce a false signal due to sparse sampling of 
the light curve. 

4. For the remaining days of overlap, the weighted differ- 
ence Am =< i«a - '«b > is determined and the goodness of 
fit estimator 
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N-2 



z 

i=l 



(m A (tj) - m B (tj + At) - Am)' 



(1) 



(corresponding to the time delay Af) is calculated. 



For the linear interpolation the errors are added in quadrature. 
The number of degrees of freedom is v = - 2 because there 
are two free paramaters: time delay and magnitude shift. The 
1 /(N - 2) factor penalizes solutions with a small number of 
overlap days. 

The best-fitting time delay can be determined by calculat- 
ing xl values for a range of time delays, and by choosing the 
time delay with the lowest xl value. In order to determine the 
associated measurement uncertainties we used 10000 bootstrap 
resamplings of the observed light curve, smoot hed by a trian- 
gular filter with a full width of 20 days (e.g., Kundic et all 
Il997l) . For each smoothed resampling the best-fitting time de- 
lay was determined for time delays between and 220 days 
(image A leading). This procedure is robust with respect to the 
size of the filter. We always interpolated image A because in 
this case the sparse sampling of our light curve between Julian 
Date-2452000=900 and 1 100 has a smaller impact on the deter- 
mined limits. We have verified that interpolating image B does 



2 Table 2 is being made available electronically at CDS. Columns 
1-4 contain the V - R colours for images A and B with errors, Column 
5 contains the Julian Date of the observations. 
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Fig. 5. Probability p for time delays based on our 
SBS 1520+530 light curve. The probability was calculated from 
10000 Monte-Carlo realizations of the light curve, and by in- 
terpolating quasar A each time. 95 per cent of the probability 
is contained in the four separate regions above the dashed line. 
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not change the answers. It would only weaken the determined 
constraints. 

The result of the Monte-Carlo resampling is shown in 
Fig. [5] For two-day bins of the time delay Af the number of 
Monte-Carlo light curves with a best-fitting time delay in the 
given range has been calculated. The probability p of each bin 
was calculated by dividing by the total number of resamplings. 
We find that there are several time delays that are consistent 
with our data; the 95 per cent confidence region consists of 
four separate regions. For these sub-regions we can calculate 
the average time delay and the standard deviation, yielding 
Af = (1 15.4 ± 2.1) days, (130.5 ± 2.9) days, (146.3 + 1.0) days 
and (198.8 + 1.1) days. The regions carry 33 per cent, 54 per 
cent, 2 per cent and 6 per cent, respectively, of the statistical 
weight. 

The sub-region with the largest statistical weight is consis- 
tent with the time delay Af = (130 + 3) days found by B02. For 
this time delay our data require an offset of Am = -0.83 mag. 
In the remainder of this paper we will use this time delay for 
our data. 

4.3. Microlensing 

In the top panel of Fig. [6] we show the quasar A and B light 
curves in one plot, where image B was shifted to the left by 
the time delay of 130 days and up by the magnitude offset of 
-0.83 mag (see Sect. 14.21 . The composite light curve has no 
large gaps. It shows that the quasar has been going through a 



Fig. 6. Top panel: /?-band light curves of the two quasar images 
A (filled circles) and B (open circles). The B light curve has 
been shifted by -130 days, as indicated by the arrow, and by 
-0.83 mag. Bottom panels: Difference light curve between the 
shifted and interpolated quasar images. This was calculated in 
both ways: by interpolating image A or by interpolating image 
B (see text). 

series of three small Am » 0.1 mag brightness variations that 
each lasted about 100 days. 

Microlensing in the lens galaxy would only affect one of 
the light paths to the quasar and could th us be detected as a 
residual light curve difference ( e.g., ISchmidt & Wambseansst 
Il998l (Wambsg anss et alll2000l) . In order to study whether mi- 
crolensing variations are present in our data, we calculated the 
difference between the two observed light curves for the time 
delay Af = 130 days found by B02 and the magnitude offset of 
Am = -0.83 mag. 

To calculate the difference light curve, quasar B was shifted 
in time and magnitude. The rest of the procedure is identical to 
the one described in Sect. 14.21 we calculated the difference by 
linearly interpolating the light curves of quasar A or B. The 
light curves were interpolated whenever the gap was less than 
40 days. No difference was calculated for larger gaps. The er- 
ror bars were added in quadrature. The resulting two difference 
light curves are plotted in the bottom two panels of Fig.|6] 

It can be taken from these plots that we do not detect a sig- 
nificant difference between the two quasar light curves. We can 
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calculate the goodness-of-fit estimator for the null hypothesis 
that there is no difference between the light curves on the N 
days for which the difference was calculated (yielding v = N-l 
degrees of freedom): 

2 1 Am\ 

x * = W^TAu^~- (2) 

i i 

This calculation yields xl — 1 -0 regardless of whether image A 
or B is interpolated, indicating that there is excellent agreement 
between the light curves. The data are formally compatible at 
a probability of 47 per cent with the null hypothesis. This xl~ 
procedure ignores a possible temporal correlation of the data, 
but there is also no evidence in Fig.|5]for such a correlation. 

Assuming that the data have a Gaussian scatter around zero 
we can determine the standard deviation 

cr 2 = — !— V Ami (3) 

1 

This yields cr = 0.03 mag if image A is interpolated and 
cr = 0.04 mag if image B is interpolated. These values are indi- 
cated in Fig. |6] (dashed lines). They are of the same magnitude 
as the error bars, again showing that the difference light curve 
is consistent with being entirely due to measurement uncertain- 
ties. 

In their earlier data taken between February 1999 and May 
2001, B02 did find a difference between the light curves of the 
two quasar images in this system. Applying our procedure to 
the light curves in their table 2, we can also calculate the differ- 
ence light curve of their data. The result is shown with in Fig.0 
(open circles) together with our difference light curve (filled 
circles) from the first of the two bottom panels in Fig.|6](image 
A was linearly interpolated). The difference curve in this figure 
is plotted without offset. 

Fig.0shows that B02 observed a coherent and highly sig- 
nificant difference light curve with a maximum amplitude of 
Am « 0.08 mag. Since the time of the observations by B02 the 
magnitude difference between the quasar images has increased 
by 0.14 + 0.03 mag (see the discussion in Sect. [5}- B02 already 
identified a linear trend (dashed line) in their data. Our data 
are consistent with this linear trend having continued until the 
epoch of our observations. 

Although the exposure times for the Maidanak and the B02 
data are similar, the error bars of the Maidanak data are larger 
than the error bars obtained by B02. The main reason for this 
are the different telescope apertures (1.5m at the AZT-22 vs 
2.5m at the NOT). In addition, however, the AZT-22 transmis- 
sion was reduced by about one magnitude before the mirror 
was cleaned on Julian Date-2452000=880 (see Fig.Q}. 

We note that the long-term variation of the difference light 
curve may even be slightly larger than shown in this plot be- 
cause some light from the lensing galaxy could be included 
in our magnitude estimate of image B (see Sect. [3}- If all of 
the lensing galaxy light were included, the Maidanak differ- 
ence light curve would have to be shifted upward by a constant 
offset of -0.08 mag to correct for the galaxy contribution. 
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Fig. 7. Composite of the difference light curves based on our 
Maydana k data (filled circle s, see Fig. |6j and the data pub- 
lished bv lBurud et all §002) (open circles, the linear trend de- 
termined by them is shown with a dashed line). In both cases 
image A was linearly interpolated. 

5. Summary and discussion 

We have presented V-band and /?-band photometry of the 
gravitational lens system SBS 1520+530 taken at Maidanak 
Observatory in the years 2003 (April to October) and 2004 
(January to August). During the « 500 day observation period 
with 80 data points in V and 123 data points in R we find small 
amplitude intrinsic variations (Am « 0.1 mag) on time scales 
of about 100 days in both quasar images. The V — R colour of 
the quasar is consistent with being constant during the observed 
period. 

Using linear interpolation of the quasar light curves, we 
have determined the 95 per cent confidence region of time 
delays for our data set. Due to gaps in the light curve, our 
data allow four separate values of the time delay, the best 
one of which agrees with the time delay of (130 + 3) days 
found by B02. Image A is leading, which is also consis- 
tent w ith lens m odels of the system (Asano l2000HFaure etall 
2002:B02 lZheleznvak et al.l2003b . 

Using the B02 time delay, we calculated the difference light 
curve between the two quasar images. This shows that within 
the statistical uncertainty the two quasar light curves are iden- 
tical during our observing interval (Fig. [6}. In the observations 
taken by B02 between 1999 and 2001, a highly significant and 
variable difference with an amplitude of Am a 0.08 mag was 
present for a fraction of the observing interval. Since then, the 
overall R-bmd magnitude difference between the A and B light 
curves has changed by 0.14 + 0.03 mag, the difference being 
larger in our data (Fig.0. 

Any variable difference between the light curves of 
SBS 1520+530 can be interpreted as gravitational microlens- 
ing because other changes of the source would be visible 
in both quasar images, delayed by the time delay. In addi- 
tion to the microlensing variability on short time-scales (« 
100 days and less) found by B02, our data show that there 
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are also variations on longer time-scales of At « (100 - 
1000) days. This overall level of microlensing variations in 
SBS 1520+530 app ears comparable to variations seen in other 
lens systems (e . g.. iHiorth et all 120021 IWvifhe & Loebll2002t 
ISchechteretaill2003l) . 

An exciting prediction for the microlensing effect of 
quasars is the colour-dependence of the microlensing light 
curve in the vicinity of caustics JWambsganss & Paczvnskil 
Il99lh . In such a situation the difference between the V and the 
R light curve could provide valuable clues to the source struc- 
ture of the quasar. We will continue to observe SBS 1520+530 
from Maidanak observatory because frequent sampling of the 
source remains crucial to derive limits on microlensing vari- 
ability. If colour variations associated with microlensing could 
be proven in this system, there would be a strong case for par- 
allel spectral observations of the quasar (see also B02). 

Since microlensing currently remains the only technique 
with the promise to scan the continuum emission regions of 
quasars on microarcsecond scales, SBS 1520+530 should be 
viewed as a prime target because of the combination of known 
quasar variability and the at least occasional occurence of mi- 
crolensing diagnostics at the same time. 
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